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Preface 


This document contains the proceedings of the 32nd annual NASA Aerospace Battery Workshop, 
hosted by the Marshall Space Flight Center on November 16-18,1 999. The workshop was attended 
by scientists and engineers from various agencies of the U.S. Government, aerospace contractors, 
and battery manufacturers, as well as international participation in like kind from a number of 
countries around the world. 

The subjects covered included nickel-hydrogen, nickel-cadmium, silver-zinc, and lithium-ion 
technologies. 
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Introduction 


The NASA Aerospace Battery Workshop is an annual event hosted by the Marshall Space Flight 
Center. The workshop is sponsored by the NASA Aerospace Flight Battery Systems Program which 
is managed out of NASA Glenn Research Center and receives support in the form of overall 
objectives, guidelines, and funding from Code S, NASA Headquarters. 

The 1999 Workshop was held on three consecutive days and was divided into five sessions. The 
first day consisted of a General Session and a Nickel-Hydrogen Session. The second day consisted 
of two Lithium-Ion Sessions with the first two papers (nickel-hydrogen and silver-zinc) placed in 
the first session due to scheduling constraints. The third and final day was a second Nickel- 
Hydrogen Session. 

On a personal note, I would like to take this opportunity to thank all of the many people that 
contributed to the organization and production of this workshop: 

The NASA Aerospace Flight Battery Systems Program, for their financial support as well 
as their input during the initial planning stages of the workshop; 

Huntsville Hilton, for doing an outstanding job in providing an ideal setting for this 
workshop and for the hospitality that was shown to all who attended; 

Marshall Space Flight Center employees, for their help in mailing the various 
correspondence, registering attendees, handling the audience microphones, and flipping 
transparencies during the workshop. 

Finally, I want to thank all of you that attended and/or prepared and delivered presentations for this 
workshop. You were the key to the success of this workshop. 


Jeff Brewer 

NASA Marshall Space Flight Center 
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^‘Standard Ni-Cd’ is a conventional Ni-Cd cell, while ‘Advanced Ni- 
Cd’ has negative electrodes with organic plate treatment. 

2 ) CPV Ni-H2 is going to be launched to demonstrate its performance 
in the orbit. 



# NASDA CPV Ni-H 2 Battery 
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# IMASOA CPV Ni ’ H 2 Battery 
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Demonstration Model 


• NASDA CPV Ni-H 2 Battery 
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Inside Structure of CPV Specification 


# NASDA Advanced Ni-Cd 
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$> NASDA Advanced Ni-Cd 
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We observed rather stable performance in the test with DOD=40%, while we found two ‘dead 
short’ cells in the test with DOD=50% at ca. 13,500 th cycle. 

When the test with DOD=40% passed ca. 15,000 cycles, we checked the V/T conditions as respect 
to DOD and temperature, (area A). We have now switched back to the test with DOD=40%. 


NASDA Advanced Ni-Cd 
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These figures show the recharge ratio (C/D ratio). The V/T level was always changed to 
keep the C/D ratio between 1.0 and 1.1. 
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0 NA 5 DA Lithium Ion Secondary Cell 
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# NASDA Lithium Ion Secondary Cell 
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• MELCO is going to establish the battery system including the charge method. For example, MELCO studies the 
electric circuits, safety unit, etc. It will also evaluate the safety of the battery system. 

• JSB checks the single cell performance. Life cycle test, safety test of cells, storage effect are included in its study. 

• NASDA focuses the evaluation on the long term performance as a battery. 


# NASO A Lithium Ion Secondary Cell 
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1) The charge method is so-called CC-CV (constant current - constant voltage) control. 


N A5DA Lithium Ion Secondary Cell 
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After ca. 170 th cycle, one of the cells started to show the lower EODV and EOCV(area A). The reason for 
these lower voltage seemed to reflect the difference of state of charge (SOC). In order to equalize SOC 
condition, we discharged each cell down to 2.75 V at 300 th cycle. After this treatment, the battery showed a 
better performance, again. 

At 500 th cycle, we decreased the ambient temperature (Ta) by 3°C, which decreased the EODV slightly. 


NASDA Lithium Ion Secondary Cell 
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This is one of the cycle data before the SOC control. A single cell showed both lower 
EOCV and EODV. 


> NASDA Lithium Ion Secondary Cell 
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This is one of the latest cycle performances of the battery with DOD=25%. After the 
treatment at 300 th cycle, the difference of voltage was reduced. 


NASDA Lithium Ion Secondary Cell 
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These figures show the temperature change through one cycle. 
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decreased the ambient temperature (Ta) by 3°C. 


# N ASDA Lithium Ion Secondary Cell 
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These results were obtained at 280 ,h cycle. In this measurement, cells connected in series showed the same 
performance. The cells kept the same SOC conditions. 


> N ASDA Lithium Ion Secondary Cell 
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These figures show the temperature change through one cycle with DOD=40%. 
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Cycle Number 

Fig. The trend data of with DOD=80%, GEO. 
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Charge and discharge curves of GEO cycle test are shown above. Stable performance could be 
observed. 


•HI NA5DA Lithium Ion Secondary Cell 
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These figures shows the temperature change through one GEO cycle with DOD=80%. 
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Above battery/cell is under development by the cooperation of Nissan Motor Co. and Shinkobe Electric 
Machinery Co. We will start our tests from spring. The point of this cell is the electrode material. 
LiMn0 2 and non-graphite carbon is used for the positive and negative electrode, respectively. 
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The cycle data of 100 Ah lithium ion secondary cell were shown. Two modes 
of LEO and one mode of GEO cycle test are to be continued. 
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capable of pulse power applications that weighs 
less than a conventional source sized for the 

same application. 


Proof of Concept Testing 
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Battery Only - 15 kLb Test 
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Ultracapacitor extends the run time of 
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Ultracapacitor enhances battery performance 

at low temperature 
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Thermal Vacuum Test Results 
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Preliminary Findings 

• Difference between methods is significant 
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Low Temperature Performance of Two Nickel Hydrogen Designs 

By 

Dr. Dean W. Maurer 1 
Loral Skynet 
PO Box 7018 
Bedminster, NJ 07921 
And 

Dr. Dawn M. See 2 

Imperial College of Science, Technology and Medicine 
Dept, of Materials 
Prince Consort Road 
London SW7 2BP, UK 


Introduction 

The Telstar 5, 6 & 7 spacecraft owned and operated by Loral Skynet and built by Space 
Systems/ Loral utilize the patented 3 Loral “Cold Charge” protocol in their battery 
operation. This protocol provides an approximately 25% capacity enhancement by 
charging at a very low temperature and subsequently discharging while ramping the 
temperature up from the end of charge temperature to as high as 30 °C. It was of interest 
to measure the relative effects of the low temperature charge and the temperature-ramped 
discharge on flight-quality cells of the two different cell designs used in the Telstar 4 
program where this protocol is not used. The Telstar 401 cell is a conventional 3.5” 
COMSAT design (50 AH) with asbestos separators and the Telstar 402R cell is a 3.5” 
modified Mantech design (45 AH) with a single layer of Zircar separator. 

The cell cycling facility used for the experimental work was designed and built some 
years ago by the AT&T Bell Labs (when Skynet was part of AT&T) for life testing of 
Telstar battery cells. When Skynet became part of Loral, this facility was donated to the 
University of South Carolina (USC). The work was carried out by one of us 2 as part of a 
graduate study program. The data were also intended for use in model verification of the 
USC first-principles nickel-hydrogen cell model. 

Experimental 

After selection of the flight battery cells, a number of the remaining cells from the 
production lots were selected for this work. The cells were placed in a temperature 
chamber in a horizontal position with a thermocouple attached to the cell midsection. 
Some of the selected cells included strain gages. The 401 cells were cycled three times at 
each test temperature. The cycle consisted of a charge for 16 hours at C/10 and a 
discharge at C/2 to 0.8 volts followed by a 3 hour open circuit stand to equilibrate the 
temperature. The 402R cells had two cycles at each temperature. First they were charged 
at C/10 for 20 hours and discharged at C/2 to 0.8 volts at constant temperature and then 

1 presenter 

2 formerly with University of South Carolina 

3 J. Hall, US Patent #5,395,706 (1995) and J. Lenhart, J. Hall and A. Applewhite, US Patent # 5,617,006 
(1997) 
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cycled again with the temperature increased during the discharge from the end-of-charge 
temperature to about 30 °C by the end of the discharge. Figure 1 is a plot of the 
temperature and pressure for the test run at -20 °C. The pressure fluctuation which was 
observed on the ramped-temperature test is an artifact caused by voltage noise on the 
strain gage wiring; this was eliminated on subsequent tests. Note that the temperature 



Fig 1 Pressure and Temperature Responses for the 402R Cells at -20 °C 

chamber was not able to maintain the temperature during the overcharge portion of the 
charge. In addition, the pressure is essentially flat after 16 hours of charge, which 
indicates minimal capacity difference between the two cycling procedures. 

Results 

The discharge curves for the test shown in Figure 1 are shown in Figure 2. The 
temperature ramping results in a “normal” discharge curve such as one would expect at 
higher temperatures, whereas the “constant” temperature discharge has the voltage 
depressed during the second half of the discharge as well as slightly less capacity. Figure 
3 shows the ramped temperature discharge curves for all of the 402 R tests. 
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Figure 2 Discharge Curves for the 402R Cells at -20 °C 

Note the uniformity of the discharge plateau voltages. The capacity achieved for 
each test is representative of the charging temperature. The results of these tests is 
further summarized in Figure 4 where the capacity is plotted vs. temperature for 
each of the temperatures in the 402R tests. The ramped discharge provides only 



Figure 3 Ramped Temperature Discharge Curves for All 402R Tests 

~3% greater capacity than the “constant” temperature, but that the difference in capacity 
between +20 and -20 °C is of the order of 50%! The difference between a more usual 
operating temperature of +10°C and the -20 °C data is still 25%. Similarly the effect 
shown in Figure 2 is also presented in Figure 5 in terms of the specific energy vs. 
temperature. In this case, there is a 5% difference at -20 °C, and a 40% energy 
density advantage between +20 and -20 °C. 
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Figure 4 Ramped vs. Constant Temperature Capacity for 402R Tests 



Temperature (°c> 


Figure 5 Specific Energy vs. Temperature for 402R Cells 


Finally the 402R cell data and the 401 cell data are compared in Figure 6. In this case the 
direct result is shown as well as a curve in which the 401 cell data are “normalized” by 
the ratio of the nameplate capacities of the cells. The only justification offered is that the 
nickel electrodes are approximately the same area in each design and the cell capacities 
differ only by the number of electrodes. 
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Figure 6 Capacity Data for 401 and 402R Cells vs. Temperature 
In Figure 6 it can be seen that: 

1) The capacity of the 401 cells below about 5 °C is relatively flat but above that 
temperature it has about the same slope as the 402 R cells. Presumably this is 
due to the decreased conductivity of the asbestos separator at lower 
temperatures. 

2) The unadjusted capacity of the asbestos cells at -20 °C is about the same as 
the capacity of the Zircar cells with fewer electrodes. 

3) The “normalized” capacity of the 401 cells at higher temperature is about 
equal to the 402R cells but the 402R cells are ~12% higher at the low 
temperature. 

Conclusions 

These studies support the generally held expectation that the Mantech cell design 
offers some performance advantages over the original COMSAT design (the COMSAT 
design is after all nearly 25 years old). More importantly the studies show a pronounced 
advantage in operating at very low temperatures: a 25% advantage at -20°C compared to 
the more usual +10°C. How this temperature control is achieved in orbit is left as an 
exercise for the thermal engineer. 
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> RR = 105%, Total Cycle Time = 90 Minutes 
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+ All cells are in pressure vessels. 
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LAUNCH AND EARLY ORBIT PHASE SIMULATION 
TESTS ON NICKEL HYDROGEN BATTERIES 

P. J. JOHNSON *, S. W. DONLEY and B. STARRITT 

TRW Space & Electronics Group 
One Space Park, Redondo Beach, CA 90278 

*E-mail: philip.johnson@physics.org 


Abstract. 

A series of experiments has been carried out with the objective of 
evaluating the performance of nickel-hydrogen batteries for a variety of 
scenarios that may be encountered pre-launch and during the early 
phases of launch and transfer orbit. Temperature of the battery was 
found to be a critical parameter. Excessive overcharge was found to 
be of no benefit for subsequent battery performance. 
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INTRODUCTION 



This presentation describes a series of experiments 
designed to elucidate the behavior of nickel- 
hydrogen batteries subject to various launch site 
and early orbit phase scenarios. The on-station 
performance has been the subject of a series of life 
tests, which have been reported elsewhere [1]. 


There were several purposes behind these test: to gather information 
that could be used to accurately determine the battery capacity 
available at liftoff (Charge Stand and Capacity Maintenance Tests), to 
determine the battery behavior and effectiveness of charging at a 
specific voltage (Transfer Orbit Operation Simulation), and to optimize 
the pre-launch battery charge operations (Launch Preparedness 
Simulation and Trickle Charge During Cool-Down). The liftoff capacity 
determination and the transfer orbit operation data are used to estimate 
the energy available at spacecraft separation to support deployment 
contingencies. 
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PHASES 



• Launch Pad 

This covers the period from the final pre-launch 
recondition to ignition 

• Ascent 

This covers the period from ignition to 
spacecraft separation from the launch vehicle 

• Transfer Orbit 

This covers spacecraft separation, solar array 
deployment, and station acquisition 


These three phases cover the time from the final on the ground charge 
to the first post eclipse charge when the spacecraft is on-station. 
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• Charge Stand Matrix 

• Launch Preparedness Simulation 

- Temperature 

- Charge return 

• Transfer Orbit Operation Simulation 

• Trickle charge during cool-down 

• Capacity Maintenance Tests 

- Normal Trickle Charge 

- Low-rate Trickle Charge 

- Very Low-rate Trickle Charge 


Normal rate for trickle charge is C/100; Low-rate trickle charge is C/330 
and very low-rate trickle charge is C/1000. The "Trickle charge during 
cool-down test was combined with the capacity maintenance tests. 
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The battery under test comprised of 23 series connected 98 Ah cells. 
The cells were provided by Eagle-Picher LLC of Joplin, MO. They were 
a dual stack back-to-back design with a single layer Zircar™ separator. 
The battery was designed, manufactured and tested by TRW Space 
and Electronics Group. 

The cells were assembled into an Engineering Model (EM) battery with 
full flight configuration. Each cell is constrained within an aluminum 
sleeve. Electrical insulation is provided by a layer of Cho-therm™ 
coated on both sides with RTV-11. The cell sleeves were then 
mounted on an aluminum baseplate. The 23 cells were arranged in an 
6x4 array, with the spare ‘slot’ used to accommodate battery 
electronics. 
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There are differences between the way cells and 
batteries are tested: 


CELLS: 

BATTERIES: 

Tested isothermally to within 3°C 

Cold plate held at constant 
temperature. Battery allowed to 
float thermally as it would 
operationally. 

Discharged at constant current, C/2 

Discharged at constant power as 
defined by system power budget. 

One hour open circuit stand between 

Four hour stand on trickle charge 

charge and discharge phases 

(C/100) between charge and 
discharge phases 


Cells undergo an industry standard acceptance test on receipt from the 
manufacturer. This allows comparison of a given lot of cells with 
previous lots of cells, which may have been procured for different 
missions with considerably different operational characteristics, for 
example GEO and LEO missions. The testing of cells closely followed 
the guidelines given in reference 2. 
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COMPARISON OF DISCHARGE AT 
CONSTANT CURRENT AND AT raw 
CONSTANT POWER 



The discharge was interrupted when the battery temperature reached 
27°C and the battery was allowed to cool on open circuit until its 
temperature had reached 18°C, when the discharge was resumed. 
Discharge was terminated when the battery voltage reached a lower 
limit of 24.15 V, which corresponded to an average of 1.05 V per cell. 
Note that this occurred much later than the maximum expected 
operational depth of discharge of 65%. 
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TEST 

CAPACITY / Ah 

NOTES 

Reference 20°C 

94-65 


72 hour stand 
at 20°C 

76-64 

81% of Reference 


The matrix consisted of one reference capacity cycle followed by one 
charge-discharge cycle with a 72 hour stand time. The stand time 
began once the battery had reached thermal equilibrium. The 
temperature was measured as the average of the cell sleeve top 
temperatures. In each cycle the charge and stand was held at 20°C 
while the discharge commenced with the battery stabilized at 4°C. The 
capacity returned after the 72-hour open-circuit stand at 81% of the 
reference test compares favorably with the result obtained during cell 
acceptance testing. 
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TEST 

TEMPERATURE 

CAPACITY 

/Ah 

VOLTAGE 
@ 1-32 h 
/ V 

NOTES 

START 

FINISH 

13°C 

27°C 

104-85 

28-713 

Tripped after 13-20 h on charge 

13°C 

18°C 

100-36 

28-892 

Tripped after 12-03 h on charge 

21 °C 

22°C 

69-39 

* 

Tripped after 8-33 h on charge 

21°C 

27°C 

92-58 

28.860 

Tripped after 12-00 h on charge 

25°C 

27°C 

74-70 

28-710 

Tripped after 9-50 h on charge 


* DISCHARGE TERMINATED EARLY 


The battery was charged to a preset end-of-charge temperature at the 
sleeve-top, followed by equilibration at 4°C (39°F) and discharge to 
24.15 V (the battery was fully discharged to measure capacity but 
voltage data are recorded after 79-2 minutes, corresponding to 65% 
DOD). 
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The above chart shows the delivered capacity at 4°C following a 
temperature limited charge. The temperature limit and the charge start 
temperature were variables. Capacity in excess of nameplate was only 
achieved at the lowest start temperature. Raising the cut off 
temperature from 18°C to 27°C led to only a marginal increase in 
delivered capacity. 
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TEST 

CAPACITY 

/Ah 

VOLTAGE 
@ 1-32 h 
/ V 

NOTES 

ID#1 140% 

101-18 

28-68 

Over Temperature after 13-2 h on 
charge 

ID#2 160% 

100-68 

28-60 

First (of three) Over Temperature after 
13-1 h on charge 

ID#3 185% 

105-61 

28-62 

First (of three) Over Temperature after 
13-2 h on charge. Cooling fault. 

ID#4 110% 

90-73 

28-99 


ID#3 185% R 

102-38 

28-65 

First (of four) Over Temperature after 
13-2 h on charge. R = Retest. 

ID#6 120% 

100.36 

28.89 

Performed during temperature tests 

ID#7 132% 

104.85 

28.71 

Performed during temperature tests 


In this series of tests, the battery was discharged to 24.15 V and the 
delivered capacity multiplied by the recharge ratio was then recharged 
at constant current. Note that only when the recharge ratio was less 
than 140% could the battery be recharged without interruption. 
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This chart shows the results of the launch pad simulation tests to 
determine optimum recharge ratio following a reconditioning. Optimum 
capacity is obtained with a recharge ratio of 130%. Higher recharge 
ratios necessitate interruption of the charge phase to prevent the 
battery overheating. This may explain the poorer capacity response at 
high recharge ratios. Excessive overcharge depresses the discharge 
voltage during the useful part of the discharge curve, l.e., during the 
maximum of 1.32 h of eclipse discharge. Depth of discharge also 
increases slightly to compensate for this voltage depression. 
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As the battery is fully charged all the charge current is converted to 
heat and the battery temperature rose quickly to its safety limit. 
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The response curve starting at 10% DOD is shown above. At the start 
of charge the charge is current limited. Little heat is generated when 
the charging efficiency is high and 6-5 Ah is returned to the battery 
before its voltage begins to throttle the charge current. As the battery 
is still cool, the current falls rapidly, and a steady state is reached after 
approximately one hour on charge, and now the battery resistance is 
17 0. The average battery temperature was below 8°C throughout the 
entirety of the charge phase of this test. 
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A similar response is achieved starting at 20% DOD, as shown in this 
figure. 24 minutes into charge the battery voltage restricts the current 
and by two hours into the test steady state is once again reached with 
the battery resistance at 17 Q. Following charge, the discharge cycle 
delivered 101-64 Ah. The average battery temperature was below 6°C 
throughout the entirety of the charge phase of this test. 
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TEST 

CAPACITY 

/Ah 

NOTES 

Standard 

11105 

Charge 14 h @ C / 10, trickle 1 A 

Capacity 


for 4 h, open 2 h. 

Capacity 

99- 02 

Trickle 300 mA for 42 min., 98 

Maintenance 1 


mA for 8h 48 min. 

Capacity 

84-53 

Trickle 300 mA for 42 min., open 

Maintenance 2 


24 h, 106 mA for 9h 27 min. 

Capacity 

87-48 

Trickle 300 mA for 42 min., 98 

Maintenance 3 


mA for 33 h 09 min. 


The capacity maintenance tests were preceded by a standard reference capacity test 
to demonstrate the battery was fully reconditioned after it had been on open circuit 
following the transfer orbit simulation tests. The capacity returned (11T05 Ah) 
compared favorably with the standard capacity value obtained in the preceding test 
(113-70 Ah). 

The first capacity maintenance cycle had no experimental stand period. Following 
charge, the trickle charge rate was 300 mA until thermal stability at 22°C was 
obtained. The period spent on trickle charge was recorded for use in the subsequent 
tests. Temperature was reduced to 4°C prior to discharge, and trickle charge at 98±5 
mA was applied during this cool down period. A capacity of 99- 02 Ah was delivered. 

Cycle 2 was the same as cycle 1 with the addition of a 24-hour stand at 22°C. This 
test delivered 84-53 Ah, suggesting that 14-49 Ah (14-8% of nameplate) were lost 
during the stand period. 

Cycle 3 repeated cycle 2 but with 98 mA trickle charge applied during the stand 
period. This test delivered 87-48 Ah, some 2-95 Ah more than cycle 2. As the trickle 
charge only added 0-098 x 24 h = 2-35 Ah, this difference cannot be attributed to the 
trickle charge alone. As the application of trickle charge necessitates the application 
of an over-potential to the battery, it seems reasonable to assume that the applied 
potential acted as a 'float' voltage and reduced the self discharge rate of the battery. 
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SUMMARY 



• A low start of charge temperature gives best charge 
acceptance 

• Excessive overcharge yields little benefit. Optimum 
recharge ratio for capacity is 130%, but for voltage is 110% 

• Battery polarization can prevent thermal runaway; fully 
charged batteries do not polarize sufficiently. 

• Trickle charge during battery cool-down from an over- 
temperature condition is unnecessary 
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Rate of charge is automatically controlled by charge determination 
and depth-of-discharge control software 



Battery Assembly 
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Pt Loading 8 mg/cm 
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Concentration 31 % 

Precharge Nickel 


Experimental Design 
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Cell / Sleeve Arrangement TRW 
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Bottom View of Cell Pack 
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Experimental Setup 
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Experimental Conditions 
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Baseline Capacities TRW 
Cell Potential vs Time 
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Time (Hours) 



Cell Capacity Comparison TRW 

(Nameplate Capacity = 160 Ah) 
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Time (Min) 


Pack Voltage vs. Cycle TRW 
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1000 2000 3000 4000 5000 6000 7000 8000 
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Discharge Voltage vs. Cycle TRW 
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Pressure vs. Cycle 
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Recharge Ratio & Pressure vs. Cycle 1 r\vv 
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Weight savings of 33% with same dimensions 
Weight savings for a 22-cell battery: 1.5+ kg 
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Cell centered in sleeve with Kapton shims at 
90° increments and potted with red RTV-566 
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No apparent 1 st mode was found in 1 0-2000 Hz range. 
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09:40:39 RNH-350 4.5 INCH CENTER-MOUNT SLEEVE 

02-Sep-1999 SN 162 COMP AXIS (Z) LEVEL 3 

Test Name: RNH-350-1_CELL-AXIS-20GRMS.008 
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09:09:10 RNH-350 4.5 INCH CENTER-MOUNT SLEEVE 

02-Sep-1999 CELL AXIS (Z) SN 143 AFTER LEVEL 3 

Sine Test Name: RNH-350JCELL.010 




09:52:08 RNH-350 4.5 INCH CENTER-MOUNT SLEEVE 

02-Sep-1999 CELL AXIS (Z) SN 162 AFTER 3RD LEVEL 

Sine Test Name: RNH-350_1CELL.014 
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similar resonant frequencies in each axis 
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Having Diameters of 2.5, 3.5, 4.5 And 5.5 
Inches. 
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147 AH CPV Cells, LEO Design (2.50 volts) 5 each 
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Capacities From 200 AH to 500 AH. 
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Each Capacity Measurement Is Made From a 16 Hour C/10 Charge and a C/2 
Discharge to Either 1 .00 Or 2.00 Volts (Depending on IPV or CPV). 
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Maximum Charge Voltage Was Approximately 3.1 Volts at 10°C. 




No Physical Damage 




for the 350 AH, 5.5 Inch Diameter Cell. 
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Voltage Rollover and Waviness 
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Comparison of Separation 
Performance for Four Variables in 
Rechargeable Silver/Zinc Cells 
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Figure 14. 
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OBJECTIVE 

Develop thermal model of Li-Ion cells in terms of heat generation, 
thermal mass, and thermal resistance. Intended for incorporation into 
battery model. 
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METHOD FOR ESTIMATING HEAT GENERATION RATE 
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Open-circuit Curves of Prismatic Li Ion cells 
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2.50 A Charge Current (C/10) 
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Heat Evolution During Discharge of Prismatic Cells 
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Cell#2 

3.58Ah 

to 3.0V 3.50Ah 
to 2.7V 3.60Ah 

3.98V 

to 3.0V 3.47Ah 
to 2.7V 3.57Ah 

to 3.0V 99.1 
to 2.7V 99.2 

Cell# 1 

3.56Ah 

to 3.0V 3.49Ah 
to 2.7V 3.57Ah 

3.98V 

to 3.0V 3. 44 Ah 
to 2.7V 3.55Ah 

to 3.0V 98.6 
to 2.7V 99.4 


Capacity In 
at 20°C 

Capacity Out 
at 20°C 

72 Hr Charge 
Retention 

Retention 

Percent 
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Trickle Charge 
Current 

0.002 Amp 

0.004 Amp 

0.002 Amp 

0.003Amp 

0.004Amp 

Cell 

SAFT 4 Ah 

SAFT 4Ah 

WGB 1.5 Ah 

WGB 1.5 Ah 

LTC 8 Ah 
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Test Plan 

Characterization: Charge Acceptance, Cell Resistance 
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Test Plan 

Simulated Leo Cycling 
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Characterization 

Charge Acceptance Results 
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Charge Voltage Limit (volts) 
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State Of Charge (ampere hours) 


Simulated LEO Cycling Results 


(/) 

tr 

TO 

6 

■O 

c 

<D 

v. 

+-» 

0 ) 

O) 

TO 


d> 

O) 

i— 

TO 

O 

</> 


o 

T3 

c 

LU 

O cb 

o 

ID 
CM 


o 

£ 


o 

Q. 

0 

c 

o 

c 

o 

jj) 

0 

O 


XL 

o 

TO 

0. 


TO 
O 

'q. o 

>» -Q 


O 

TO 

Q. 


— _ 0 


0 

O 


I 


0 U) 

c 


0 

1 

CM 

I 


0 

I 


</) 

t 

TO 


■D 

C 

0 


0 
O) 
TO 

O 

> 

0 

Ui 

i_ 

TO 

O 
</> 

£ 5 

0 
O 


O 

■o 
c 
LU 

O CD 

o 

ID 


O 

Q. 

0 
C 

o 

c 

o 

O "0 

To 0 

o JZ 

Q. O 

>> -Q 


O 

£ 


X 

O 

TO 

0 . 

0 

o 


XL 

O 

TO 

Q. 


0 

0 

O 

0 

U) 


0 
o 

cm '55 


1 1 1 


tfi 

</5 

>* 

TO 

c 

TO 


0 

0 

Q. 

(/> 

■ MM 

T3 

O 

TO 

Q. 

0 

O 

■ 

CD 


1999 NASA Aerospace Battery Workshop 


- 359 - 


Lithium-Ion Session II 


Wksp9901.ppt 


25% DOD LEO Cycling at 25 Deg C -- 6-Cell Pack 
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25% DOD LEO Cycling at 25 Deg C - 2-Cell Pack 
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25% DOD LEO Cycling at 25 Deg C - Single Cells 
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25% DOD LEO Cycling at 15 Deg C -- 6-Cell Pac 
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25% DOD LEO Cycling at 15 Deg C -- 2-Cell Pack 
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25% DOD LEO Cycling at 15 Deg C - Single Cells 
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6-Cell Pack EODV Dispersion Analysis 
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Summary 
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balancing) is acceptable and is not increasing. 
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Li-Ion FOR SPACE APPLICATIONS 
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SAFT Development Context 
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Prototype Cell Development 
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Prototype Cell Development (Cont’d) 
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▼ Mains evolutions from Stentor to industrial design 
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□ Electrode loading in order to improve cell cyclability 



Industrial Line 
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Li-Ion INDUSTRIAL SPACE CELL 
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80 % DOD on STENTOR Cell 
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100 % DOD Industrial Cell 
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Negative capacity loss impact 
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Computational Fluid Dynamics Modeling of Nickel Hydrogen Batteries 


R. Cullion, W.B. Gu, C.Y. Wang 
GATE Center for Advanced Energy Storage 
Department of Mechanical and Nuclear Engineering 
Pennsylvania State University 
University Park, PA 16802 

P. Timmerman 
Jet Propulsion Laboratory 
Pasadena, CA 91109 


Abstract 

An electrochemical Ni-PU battery model has been expanded to include thermal effects. A 
thermal energy conservation equation was derived from first principles. An 
electrochemical and thermal coupled model was created by the addition of this equation 
to an existing multiphase, electrochemical model. Charging at various rates was 
investigated and the results validated against experimental data. Reaction currents, 
pressure changes, temperature profiles, and concentration variations within the cell are 
predicted numerically and compared with available data and theory. 

Introduction 

The nickel hydrogen battery provides an important form of energy storage. This battery 
has applications in both the aerospace and the electric vehicle industries. Nickel 
hydrogen batteries are used in these fields due to their ability to provide long cycle 
lifetimes, high specific energy, high power density, and high tolerance to overcharge. 

Demands by industry have made mathematical modeling of these batteries crucial, in 
particular, thermal modeling has become essential. Temperature variations within the 
battery influence the performance, life, and reliability of the power source. At low 
temperatures there are capacity losses due to high internal resistance within the battery 
and at high temperatures the capacity is compromised by fast self-discharge. Use of the 
active material must be balanced in order to maximize battery life. To do this a uniform 
temperature profile is necessary. Modeling provides a way to monitor the temperature 
profiles, and a way to test improvements without running lengthy and expensive 
experiments. Another very important reason thermal modeling is essential is safety. Side 
reactions occurring within the battery often result in temperature increases, which could 
lead to thermal runaway. Modeling would allow the runaway to be predicted. 

The cell electrochemical behavior is described by electrochemical reactions that are 
affected by temperature-dependent electrochemical and transport properties. This 
interrelation of the electrochemical reactions and temperature-dependent properties 
makes it necessary to couple the electrochemical and thermal models of the NLH 2 cell 
into a single comprehensive model. Incorporating the thermal energy conservation 
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equation, in which the local heat generation rate is determined by the electrochemical 
processes, into the model achieves the coupling. 

A one-dimensional, thermal and electrochemical model has been developed by DeVidt et. 
A1 [2]. The transport phenomena of the gaseous species and electrolyte within the cell 
were described using concentrated solution theory and the volume averaging technique. 

A pseudo-second dimension was used to model the proton diffusion inside the active 
material. Three reactions in the nickel electrode were considered. The main reaction was 
the oxidation/reduction of the nickel species into the solid active material, and the side 
reactions were oxygen evolution/reduction and oxidation of dissolved hydrogen. The 
present work provides a two-dimensional, fully coupled, thermal and electrochemical 
model. It also makes use of the volume averaging technique and concentrated solution 
theory [8]. Unlike the model developed by DeVidt et. al. [2], the present work employs a 
multiphase reaction scheme. At the nickel electrode two main reactions are considered. 
Oxidation/reduction reactions of the nickel species in both the beta and gamma phases 
are considered. The side reactions, oxygen evolution/reduction and oxidation of dissolved 
hydrogen, are also included. The gamma phase is less stable than the beta phase and 
occurs mainly during overcharge [6]. Adding both the beta and gamma phase reactions to 
the model allows for more accurate prediction of the discharge/charge behavior as was 
demonstrated by Timmerman [7]. 

DeVidt’s work [2] considers the cell to be isothermal spatially. An adiabatic case and an 
isothermal case were investigated. The focus of this study is to create an axisymmetric, 
multidimensional, thermal and electrochemical coupled model. The multidimensional 
nature of the present model is necessary because temperature gradients in the Ni-Fh cell 
are present mainly in the radial direction due to heat dissipation through the battery wall, 
while the electrochemical and ion transport processes occur in the axial direction. 
Temperature gradients within the cell will be determined by this model due to their 
importance in predicting cell performance, making improvements, and for safety 
purposes. The energy equation of the two-dimensional model will include a heat 
generation term that incorporates the heat effects due to electrochemical reactions, joule 
heating, and phase changes of non-electrochemical nature. Also, as mentioned earlier, 
DeVidt’s model [2] considered two extreme boundary conditions, an adiabatic case and 
an isothermal case. This model allows for a general convective boundary condition. 

Thermal-Electrochemical Model 

A single Ni-H2 cell consists of an aluminum current collector, a nickel electrode, a 
separator, a hydrogen electrode (which has a micro-porous Teflon membrane backing 
permeable only to gas), a nickel current collector, and a gas diffusion screen, as shown 
schematically in Fig. 1 . A concentrated KOH aqueous solution serves as the electrolyte. 
The Ni-H 2 battery is composed of a stack of these disk-shaped cells housed within a 
pressure vessel [3]. The cells are center aligned within the vessel. There is an outer gap 
between the wall of the pressure vessel and the circumference of the cell stack to ensure 
electric insulation between the cell stack and the pressure vessel and to provide a path for 
hydrogen gas to diffuse from the external reservoir into the cell stack. It is important to 
note the cell is axisymmetric, therefore although only the radial and axial directions will 
be examined the results will be the same as those found for a three-dimensional case. 
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Figure 1. Schematic Representation of a Ni-FF Cell 


Assumptions 

The present work is based on the following assumptions. This problem is considered to 
be axisymmetric. The convection effect between the outer edge of the cell and the 
pressure vessel wall, important when the temperature difference between the two is large, 
is ignored for simplicity. As a result, in this space the species transport is by diffusion and 
migration only and the heat transfer is by conduction only. The nickel electrode consists 
of composite cylindrical needles with a substrate inside (see Fig. 1). The active material 
film is assumed to be a mixture of multiple phases with a constant diffusion coefficient 
for the protons. Furthermore, the a and y phases are treated as a single phase [7]. The 
electrode porosities are considered to be constant. There is a continuous gas-pore network 
throughout the cell with a uniform constant volume fraction. The gas, a mixture of 
hydrogen, oxygen, and water vapor, is ideal. Its composition varies with time but does 
not change with position [1], The solid phase is assumed to be completely wetted by the 
electrolyte film. Interfacial chemical equilibrium exists in the liquid phase for all species 
other than the dissolved oxygen that has a small value of mass diffusivity in the liquid 
electrolyte. Electrical equilibrium exists in the liquid phase due to the large value of ionic 
conductivity of the electrolyte. Finally it can be said that the combined thickness of the 
current collector and the Ni electrode is much smaller than the diameter of the cell, 
therefore <f> SU b is considered to be constant along x but to vary with r. 

Model Equations 

At the nickel electrode the following main reactions occur: 

(5 phase reaction: 
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mOOH p + H 2 0 + e~ <» Ni(OH) 2 „ + OH~ 
y phase reaction: 

NiOOH r + H 2 0 + ne~ <=> Ni(OH) 2 a + OH ' 
y phase production: 


(la) 


(lb) 


NiOOHp + ne~ + K + -> NiOOH r + O/T 


(lc) 


a phase reconversion: 

M(0//) 2 . a -> M(Oi50 v 


(Id) 


Oxygen evolution/reduction and the oxidation of dissolved hydrogen also occur at the 
nickel electrode and are considered to be side reactions: 


-0 2 +H 2 0 + 2e- <=>20H- 
2 2 2 


2H 2 0 + 2e~ <— H 2 + 20//' 


( 2 ) 

(3) 


Two reactions occur at the hydrogen electrode, the main reaction: 


H 2 +20H- o2// 2 0 + 2e‘ 


(4) 


and the side reaction is: 


2 OH' <--0 2 +H 2 0 + 2e 
2 2 2 


(5) 


The general Butler-Volmer equation is used to determine the rates at which these 
reactions occur. The equations are written with respect to a specified reference state, 
namely a temperature of 10°C and an electrolyte concentration of 7. 1 d-3 mol/cm 3 . 
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The subscripts in the above equations; la, lb, lc, 2, 3, 4, 5 refer to the reactions labeled 
with the respective equation number. 


The overpotential of a given reaction j, qj, is given by [4]: 

( ( \[ dU dU Kf \\ 

'l J = K-t-\U J , re/ + (T-Tj—- 


cT dT 


(13) 


The temperature dependency of the overpotential is considered in the above equation. 
The reference overpotentials for the nickel electrode are calculated using the following 
equations [7]: 


„ RT 

U^ nf = U:.-— lnl 


^ sref 

- r H 

\ u 5max sref y 


+ f^-l) 


Uu* = K 


(14) 

(15) 


and the reference overpotentials for the oxygen and hydrogen reactions are [1]: 


RT ( \ 

U>,ef = U 4j ref = - 0.9263 — ln(/7"; ) (16) 

u 2 „, = l/ w = 0.3027 + — \r{p % ) (1 7) 


k in equation 14 is an intercalation constant used to describe the interaction between 
species in a solid solution. 


The conservation of OH species in the liquid phase is represented by the following 
equation [7]: 



dt 


v(d""Vc°")+ 


t°_ -1 

F 


• OH 


(18) 
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where t°_ is the transference number of OH' with respect to the velocity of the solvent. 
j OH is the total transfer current from all electrochemical reactions that generate or 
consume OH at the electrode/electrolyte interface: 


f. OH .OH 
J b +J g 


in the nickel electrode 


j 0H = S 


0 


in the separator 


< aH(T n 4 + Tn5) in the hydrogen electrode 


(19) 


where 


j OH p- aNi(l - X y )[T nla + 0 n2 + t n3)p] (20) 


j = a Ni X y [I n 1 b + (l n2 + T n3) y ] 

Similarly, j 2 and j 2 are the transfer currents associated with the generation or 
consumption of hydrogen and oxygen respectively: 


a NHn3 



in the nickel electrode 
in the separator 


< a HHi4 in the hydrogen electrode 


( 21 ) 


aNUn2 



in the nickel electrode 
in the separator 


^a{-[Tn5 in the hydrogen electrode 


( 22 ) 


Conservation of H 2 species in the liquid phase is given by [7]: 

d ( £eC " 2 ^ = V(£# Vcf 2 )+ 2- j H > + J" g 


dt 


IF ' 


O 2 species conservation in the liquid phase is [7]: 

= v(dS Vcf- )+ ±j°' + J% 


4F 


(23) 


(24) 
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The interfacial mass transfer rate of the hydrogen or oxygen from the liquid to the gas 
phase is expressed as: 


Hi 


RT C g ' C ej 


i = H 2 and 02 


J = K» 

eg V 

Conservation of species in the gas phase is given as [7]: 




dt 


Conservation of species in the solid phase [7]: 


akcf) j ; 

dt F 


( 25 ) 


(26) 


(27) 


where 


aNi(l -Xg)T n la for5 = b 

aNiXglnlb for5 = g 

Conservation of charge in the liquid phase is represented by [7]: 

V ■(k eff V<t > e )+ v(kf V In c OH )+ j OH =0 
Conservation of charge in the solid phase is given [7]: 



(28) 


(29) 


v(<rf ■ V<0-yf" = 0 

”sb 

The energy conservation equation is written in the following form [4]: 

pC D — = V / NT + q m 
8t * 


where 


PC P = £ ePe C p + £ sPs C p + £ gPg C p 

A = s e A e f + s s A f + £ g A*g 

The third term of the energy equation has three components. The first part: 


(30) 


(31) 


(31a) 

(31b) 
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(32) 


pKZ',fc+ n>r 

represents the heat effect produced by the electrochemical reactions [4]. This component 
includes both the reversible heat characterized by the Peltier coefficient: 

dU , 

n ,=T — J - (33) 

J dt 

and the irreversible heat related to the surface overpotential. 

The second component of the heat generation term is associated with phase changes of 
the non-electrochemical nature, e.g. evaporation or condensation of water within the cell 
[4]: 


£A/,r = AVv + ^<"‘ 0 - / ."' < ’ 

The third component of the heat generation term: 


(34) 


7T +kfvinc 0l < Vt,)dV (35) 

* y a 

represents the joule heating in both the solid and liquid phases [4], 

Electrochemical and transport properties that are temperature dependent are expected to 
follow the Arrhenius equation [4]: 


O = <t> 0 exp 



f 

^ act 

f 1 

Ml 


R 

[t 

~T 0 )) 


(36) 


where E act is the activation energy corresponding to a property O. 

The top and bottom of the cell were considered to be adiabatic. The thermal boundary 
condition applied to the inner edge of the cell and the outer wall of the pressure vessel 
was the following: 

~ k % = h(<T ~ ) ( 3? ) 
where k is the thermal conductivity of the material, and h is the heat transfer coefficient. 

Results and Discussion 

Cell performance under various conditions was simulated using the model described 
above. Physical properties and electrochemical parameters were provided from 
experimental work conducted at JPL and are listed in Table I. Several parameters were 
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not available, but found through trial and error to produce the best overall agreement with 
experimental charge data (these values are indicated with an asterisk). 


Table 1: Parameters 


Parameter 

Values 

References 

Kinetics 

io,la 

3.0d-04 A/cm 2 

* 

io,lb 

1 .0d-04 A/cm 2 


io,ic 

2.0d-55 A/cm 2 


io,2 

3.0d-10 A/cm 2 

* 

io,3 

2.5d-16 A/cm 2 

* 

io,4 

6.5d-05 A/cm 2 

* 

io,5 

1 ,0d-07 A/cm 2 

* 

aia 

5.0d-01 A/cm 2 


aib 

5.0d-01 A/cm 2 


au 

5.0d-01 A/cm 2 



1 ,0d+00 A/cm 2 


a 3 

5.0d-01 A/cm 2 


a4 

5.0d-01 A/cm 2 


0.5 

1 ,0d+00 A/cm 2 


P hydrogen, ref 

52.0 Atm 

Ref 1 

P oxygen, ref 

5.0e-4 Atm 

Ref 1 

U, a 

4.2d-01 V 


u Ib 

3.35d-l V 


U, c 

4.4d-01 V 


u 2 

3.0d-01 V 


U 3 

-9.26d-l V 


u 4 

-9.26d-l V 


u 5 

3.0d-01 V 


Ela 

5.0d+04 J/mol 


Eib 

3.0d+04 J/mol 


Eic 

0.0d+00 J/mol 


e 2 

8.0d+04 J/mol 


e 3 

0.0d+00 J/mol 


e 4 

2.0d+04 J/mol 


e 5 

1 .2d+05 J/mol 


Thermal 

dudtu 

-2.1 Id- 3 V/K 


dudtib 

-1.51 d-3 V/K 


dudtic 

0.0d+00 V/K 


dudt 2 

-1.68d-3 V/K 


dudt 3 

-8.36d-3 V/K 


dudt* 

-8.36d-3 V/K 


dudts 

-1.68d-3 V/K 


initial temperature 

15°C 


reference temperature 

15°C 


environmental temperature 

0°C 
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heat transfer coefficient 

5 W/m 2 K 

Cell 


electrode pairs 

17 

inner pressure vessel radius 

4.435 cm 

outer electrode radius 

4.400 cm 

inner electrode radius 

0.350 cm 

A1 current collector thickness 

0.0015 cm 

Ni electrode thickness 

0.453 cm 

Separator thickness 

0.210 cm 

Pt H 2 electrode thickness 

0.400 cm 

Ni current collector thickness 

0.0010 cm 

Gas screen thickness 

0.0610 cm 

Porosity of substrate 

0.85 

Porosity of Ni electrode 

0.41 

Porosity of Pt H 2 electrode 

0.85 

Porosity of separator 

0.85 

Specific interfacial area of 


substrate 

2000 cm 2 /cm 3 

Specific interfacial area of 


nickel electrode 

3919 cm 2 /cm 3 

Specific interfacial area of 


Pt-H 2 electrode 

2000 cm 2 /cm 3 

H 2 reference concentration 

1.461d-3 mol/cm 

O 2 reference concentration 

1.529d-3 mol/cm 

Reference concentration 

7. 1 d-3 mol/cm 3 


Experimental data for charging at various rates and a constant temperature of 1 0°C were 
used to validate the present model. Figure 2 compares the experimental cell potential 
curves with those predicted. The general shape compared well. The quick rise during 
initial charging was captured, and the second plateau was reproduced. At lower rates 
there is an under-prediction of the potential and as the rate increases there is a slight 
overprediction. Because a single set of parameters was kept for all simulations the 
observed discrepancies were expected. 

Figure 3 compares cell pressure variations measured experimentally with the present 
model predictions. It plots the cell’s pressure against time. Excellent agreement is shown 
for all three charging rates. The change in cell pressure is approximately linear until the 
overcharge period begins. At that point the pressure rise becomes strongly non-linear. 


1999 NASA Aerospace Battery Workshop 


- 526 - 


Nickel- Hydrogen Session II 


.7 



Fig. 2: Comparison of predicted and measured charge curves for a Ni-Hh cell at 
different charge rates and constant 10°C. The lines are simulation results, and the 
symbols denote experimental data. 



Figure 3. Comparison of experimental and simulated pressure variations for 
charging at various rates at 10°C. 
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Reaction current [A/cm3| 


The following results obtained from the computer simulation are intended to further 
reveal the internal processes of the Ni-H 2 cell during charge. 

Figure 4 plots the reaction current distributions within the nickel electrode when the cell 
is being charged at a rate of C/1 and a constant temperature of 10°C to 50% charge input. 
As discussed earlier there are both primary reactions and secondary reactions that occur 
in the nickel electrode. Theoretically, prior to reaching 100% charge input the primary 
reactions are the principle reactions. Beyond 1 00% charge input the secondary reactions 
take over as the major reactions. Figure 4 clearly shows that for 50% state of charge the 
primary reactions, la and lb, are dominant. It is also important to note that these 
reactions are stronger at the electrode/electrolyte interface than at the back of the 
electrode. Inspection of the chemical reactions, la and lb, reveals that OH' is a reactant 
during charging. A large supply of OH' explains why the reaction is more intense at this 
location. 



Figure 4: Reaction current distribution inside the nickel electrode for a cell charged 
at C/1 at constant 10°C to 50 % charge input. 

Figure 5 also shows the reaction current distribution within the nickel electrode when the 
cell is being charged at C/1 and kept at 10°C. It is charged to 150% charge input in this 
case. As expected, the simulation shows that the primary reactions diminish to 
approximately zero and the secondary reactions develop once the battery has been 
charged to 100% charge input. The oxygen and hydrogen reactions, 2 and 3, become the 
dominant reactions once the 100% point is passed. Also, from Fig. 5 it can be noted that 
the secondary reaction rates are larger at the electrode/electrolyte for the same reason as 
stated previously. 
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Figure 5: Reaction current distribution within the nickel electrode of a cell charged 
at C/1 at constant 10°C to 150 % charge input. 

Figure 6 plots the reaction currents in the nickel electrode as a function of the charge 
input relative to the cell capacity. The primary reactions are seen to dominate until 100% 
charge input is reached, beyond that point the secondary reactions take over. Also to be 
noted is that the total current prior to and after 1 00% charge input remains the same. 
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Figure 6. Reaction currents within a cell charged at C/1 and constant 10°C as a 
function of percent charge input. 
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Figure 7 plots the partial pressures of the hydrogen and oxygen within a cell charged at 
C/1 and 10°C. As expected, the increase of pressure due to the hydrogen reaction is 
dominant. The primary reaction at the hydrogen electrode is the main contributor prior to 
1 00 % charge input, and beyond 1 00 % the secondary hydrogen reaction at the nickel 
electrode is adding to the pressure increase. Also at 100 % charge input and beyond the 
secondary oxygen reaction at the nickel electrode begins, producing a rise in pressure and 
contributing to the overall increase in pressure. 



Time [hr] 


Figure 7: Partial pressures within a cell charged at C/1 and 10°C to 150% of the 
cell capacity as a function of time. 

Figure 8 explores the thermal effects within the cell, showing the temperature profile 
within the cell with respect to the cell radius. As mentioned earlier the cell is shaped like 
a disk with a hole at its center. The cell is hottest near the inner edge of the hole and the 
temperature decreases as the outer diameter of the cell is approached, this is expected 
based on the boundary conditions. A temperature change of 3°C is seen. A temperature 
difference of this magnitude influences the kinetic rates of the reactions by as much as 
10 %. 
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Figure 8: Average temperature of the cell charged at C/1 to 150% of its capacity as 
a function of the radius with h = 5 W/m 2 K. 



Figure 9: Average temperature of a cell charged at C/1 to 150% of its capacity as a 
function of time, when h = 5 W/m 2 K and Tj„j t j a i = 15°C. 
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Figure 9 is a graph of the change in average cell temperature with time for a cell charged 
at C/1. This plot shows an initial temperature drop during the early stages and then an 
increase in temperature during the later times. At the end of the charging period the 
battery is hotter than it was initially. The primary reactions dominate when the charging 
time is less than one hour. It is known that the primary reactions are endothermic [5], 
therefore explaining the temperature drop during the beginning stages of charging. The 
hydrogen reaction at the hydrogen electrode is exothermic as is the hydrogen side 
reaction at the nickel electrode. The hydrogen reactions coupled with resistive heating 
cause the temperature rise as charging progresses [2]. 

Experimental temperature measurements were made at the pressure vessel seam. This 
data was then plotted versus the charging time, as shown in Figure 10. An initial drop in 
temperature is followed by a strong increase in temperature. Qualitatively Figure 9, the 
simulation results for the average temperature of the cell, demonstrates a trend very 
similar to the temperature profile found through external measurements as seen in Figure 
10 . 



Figure 10: Experimental pressure vessel temperatures plotted as a function of the 
charging time for a cell charged at C/1. 

Figure 1 1 addresses two important points. It plots the cell proton concentration contour 
for a cell charged at C/1 to 50 % of the cell capacity. It demonstrates the necessity of 
having two-dimensionality in the model. The plot shows that the proton concentration 
varies in both the axial and radial directions. It also shows the effect of temperature on 
the cell. The proton depletion is greater at the inner diameter of the cell, this is related to 
the higher temperatures at this location, clearly demonstrating the thermal, 
electrochemical link. Also, as was noted earlier, there is larger proton concentration 
depletion along the electrolyte interface due to the plentitude of OH'. 
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Figure 11: Proton concentration in a cell charged at C/1 to 50 % the cell capacity. 


Conclusions: 

A two-dimensional, fully coupled, thermal and electrochemical Ni-H2 model has been 
developed. Comparison with experimental data validates the model’s results. This model 
complements experimentation by allowing the user to understand occurrences within the 
cell for various conditions and for many different charging rates. Examples of the 
program capabilities have been given. Work is underway to expand this model to the 
entire cell stack, and the thermal interactions between the cells that arise due to this 
change will be addressed and investigated. 
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List of Symbols 

a specific interfacial area, cm 2 /cm 3 

c' volume-averaging concentration of species i over a phase, mol/cm 3 

C p specific heat, J/g/K 

D' diffusion coefficient of species i in a phase, cm 2 /s 

E ac t activation energy, J/mol 

F Faraday’s constant, 96,487 C/m ol 

h heat transfer coefficient, W/m 2 K 

i n j area-averaged transfer current density of reaction j. A/cm 2 

ipj.ref exchange current density of reaction j at reference conditions, A/cm 2 
j' reaction current density due to production or consumption of species i, 

A/cm 2 

k thermal conductivity, W/mK 

q’” heat generation term 

R universal gas constant, 8.3143 J/mol K 

R S b electrical resistance from solid/substrate interface to the bulk of solid, 

Q cm 2 
t time, s 

t°. transference number of OH' with respect to the solvent velocity 

T absolute temperature of the cell system, K 

Uj, re f open-circuit potential for reaction j at reference conditions measured with 
respect to a Hg/HgO reference electrode, V 
Vj volume of a phase, cm 3 

Greek Symbols 

a a j, a C j anodic and cathodic transfer coefficients for reaction j 
Ah enthalpy change 

£i volume fraction of a phase in the reference volume fraction 

r|j surface overpotential of electrode reaction, j, V 

<f> general symbol representing a property at a temperature T 

X. thermal conductivity, W/m 2 K 

n Peltier coefficient 

<|> potential in a phase, V 

p density of a species I, g/cm 3 

a thermal conductivity W/m 2 K 


1999 NASA Aerospace Battery Workshop 


- 535 - 


Nickel-Hydrogen Session II 


Pane Intentionally Left Blank 


VALIDATION OF THE STATIC MODEL FOR 
NICKEL-HYDROGEN CELLS 


By: 

Lawrence H. Thaller 
The Aerospace Corporation 
El Segundo, CA 

Presented at: 

The 1999 NASA Aerospace Battery Workshop 
Huntsville, AL 
16-18 November 1999 


Electronics Technology Center BES THE AEROSPACE 

Energy Technology Department KsSlCORPORATION 


1999 NASA Aerospace Battery Workshop 


- 537 - 


Nickel-Hydrogen Session II 


OUTLINE 


• BASICS OF STATIC MODELS 

• MODEL FOR NICKEL-HYDROGEN CELLS 

• APPLICATIONS OF A VALID MODEL 

- Assess projected cycle life 

- Trades single layer vs. double layer separators 

• EARLIER VALIDATION EFFORTS 

• ADDITIONAL FACTORS CONSIDERED HERE 

- Density variations with DOD 

- Compressibility of components 

- Survey of available dryout information 

- A closer review of earlier partitioning studies 

• CURRENT STATUS OF MODEL 
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BASICS OF STATIC MODELS 


• A SERIES OF CALCULATIONS TO DETERMINE POROUS 
WETTABLE VOLUME WITHIN THE CELL 

- Densities, porosities, dimensions 

• INCORPORATES CHANGES DURING CYCLING 

- Cycle to cycle, start to finish 

- Electrode expansion, plaque corrosion, gamma phase formation 

• DETERMINE OR PROJECT IMPACT OF CHANGES IN 
ELECTROLYTE PARTITIONING BETWEEN COMPONENTS 


Electronics Technology Center yaffil THE AEROSPACE 

Energy Technology Department K&i CORPORATION 


A static model is one that calculates the total wettable volume of the 
components within a cell. By knowing the interior volume of a cell and 
the amount of electrolyte in the cell, an assessment can be made as to 
whether there is too much or too little electrolyte. The desired amount of 
electrolyte depends on whether the cell design is a flooded one or a 
starved one. It also depends on the changes that are expected to take 
place during the course of expected usage or cycle life. We have 
developed static models for nickel-cadmium, nickel-hydrogen, silver-zinc, 
and lithium-thionyl chloride cells. For the nickel-hydrogen model which 
will be addressed here, the model was set up to investigate the dryout 
tendencies that are caused by the natural degradation mechanisms as 
the cells are cycled. Expansion of the nickel electrode, corrosion of the 
sintered nickel plaque material, and the formation of the gamma phase of 
charged active material all result in reductions in the volume of 
electrolyte. Since the rates of these processes are partially known, 
information from a validated static model can be helpful in evaluating the 
projected cycle life of a particular cell design. 
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STATIC NICKEL-HYDROGEN MODEL 


• FULL DESCRIPTION GIVEN AT THE 1997 IECEC MEETING 
- Details of how to construct a model for your own use 


• PAPER GIVEN AT 1998 IECEC MEETING DESCRIBED 
DATABASE OF DIFFERENT DEGRADATION MECHANISMS 

• 1999 PRESENTATION AT THE ANNUAL BATTERY 
CONFERENCE DESCRIBED VOLUME TOLERANCE OF CELL 

• SOME FINE TUNING EFFORTS WERE STILL NEEDED 


Electronics Technology Center ISjl THE AEROSPACE 

Energy Technology Department ISsScOR POR ATI ON 


A full description of the initial version of the model was presented at the 
1997 IECEC meeting in Hawaii. The model consists of a spreadsheet 
wherein all the diameters, thicknesses, porosities, and densities are used 
to calculate the wettable volumes of all the components. There are only 
a few design-specific variables that are required to be inserted into the 
spreadsheet. These include the amount of plate expansion that is 
expected, the amount of plaque corrosion that might take place, and the 
amount of charged active material that will be found in the gamma phase. 
Of course, the amount of electrolyte added to the cell during activation is 
also needed. The resulting calculations will determine the amount of 
electrolyte that will be in each component given the expected changes 
that were inserted into the model. 

Extensive review of life cycle testing data coupled with post test 
destructive physical analyses have permitted reasonable estimates to be 
made of the rates of plaque corrosion and electrode expansion. A 
summary of this ongoing study as of 1998 was presented at the 1998 
IECEC meeting in Colorado Springs, CO. 

More recently, a technique was developed and reported at the 1999 
Annual Battery Conference held in Long Beach, CA. This technique 
permitted the increase in internal resistance of a nickel-hydrogen cell to 
be monitored as the volume of electrolyte was reduced to the point where 
the cell performance was unacceptably low. 

Since then, other efforts have been explored in an attempt to match the 
results of post test analyses with results predicted by the model. 
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APPLICATIONS FOR A VALIDATED MODEL 


• DRYOUT RESULTS IN AN INCREASE IN INTERNAL 
RESISTANCE 

- Gradual at first, becoming significant if diffusional limit is 
exceeded 

• DRYOUT OCCURS DURING CYCLING 

- Aggravated by aggressive cycling conditions 

• THEREFORE, MODEL HELPFUL IN ASSESSING: 

- Cycling conditions when long cycle life is required 

- Evaluation of cell designs with single layer of separator 

- Evaluation of cells that are under filled with respect to 
electrolyte 


Electronics Technology Center PgS THE AEROSPACE 

Energy Technology Department Csl C O R P O R AT I O N 


Models of any type are developed because it is hoped that results, 
information, or insight can be gained that would not be obtained using 
other means. As with any other cell type that is classified as being a 
“starved” design, nickel-hydrogen cells have been known to display 
symptoms of having an insufficient amount of electrolyte. This situation is 
usually referred to as a dryout condition. Besides an increase in internal 
cell impedance, dryout can result in severe cases of a condition that 
results in a diffusion-limited current. In this situation, there are not 
enough ions to conduct the current through the electrolyte and the result 
is a sudden drop in current carrying capability. 

Since the degradation mechanisms in nickel-hydrogen cells are 
reasonably well understood as well as the factors of temperature and 
recharge protocols that can exacerbate the rates of degradation, a 
validated model would be useful for commenting on suggested cell 
designs, cycling conditions, and recharge protocols. 


1999 NASA Aerospace Battery Workshop 


- 541 - 


Nickel-Hydrogen Session II 


MAJOR ASSUMPTIONS OF MODEL 


• ELECTROLYTE DISTRIBUTES ITSELF ACCORDING TO THE 
PORE SIZE CHARACTER OF THE COMPONENTS 
- Capillary pressures are equalized 


• THE SEPARATOR GIVES UP ELECTROLYTE TO THE 
EXPANDING NICKEL ELECTRODE 

• THE EFFECTIVE DENSITIES OF THE DIFFERENT FORMS OF 
THE ACTIVE MATERIALS ARE AS FOLLOWS: 

- Discharged beta material - 3.97 

- Charged beta material - 4.56 

- Gamma material - 3.71 


Electronics Technology Center pga THE AEROSPACE 

Energy Technology Department KS C O R P O R AT I O N 


Like most models, there were a few assumptions in the original model. 
The major one in the initial version of the model was that the wettable 
pores of the nickel electrode would remain filled witlvelectrolyte at the 
expense of the electrolyte contained in the separator. 

When two porous materials are in contact and have some amount of 
liquid within their pore structures, the liquid will distribute itself according 
to the wettability and the pore size character of the two components. 

This redistribution will result in equal capillary pressures between the 
pores of the two structures. Based on an extensive study* that was 
carried out earlier when the topic of pore size engineering was being 
developed, the pore size distribution of typical cell components were 
determined. For the cases where zircar was the separator, the difference 
in pore size distribution vs. a typical nickel electrode suggested that the 
nickel electrode would remain filled. 

The second assumption deals with the densities of the different forms of 
active material within the nickel electrode. The ones used in this study 
are slightly different than those found in the normal handbooks. They 
have been derived from practical experience with nickel elecrodes. 

The third assumption was related to the lack of compression that the cell 
components experienced while in the cell. 

In this study, the first and the third assumption will be reviewed and 
incorporated into an upgraded version of the model. 

*Abbey, K. and Briton, D. “Electrolyte Management in Porous Battery 
Components”, NASA TM 83073 
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EARLIER NASA STUDIES ON DISTRIBUTION 
AND PARTITIONING OF ELECTROLYTE 


NASA TM 83073 by Abbey and Britton - Electrolyte Management in Porous 
Battery Components 
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call components. 

<» Percent saturation of Individ- 
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The above two figures that appeared in the referenced NASA publication. 
They were generated starting with the pore size distributions of the 
individual components of a typical nickel-hydrogen cell - single layer 
zircar, standard nickel electrode, and the normal hydrogen electrode. 

With the given pore size distribution data and the requirement to have 
equal capillary pressures across the three components, these two curves 
were developed. 

The figure on the left charts how the electrolyte is partitioned across the 
three components as electrolyte is withdrawn from the set of parts. 

When the components are filled with electrolyte, about 45% of it is in the 
nickel electrode, another 45% in the separator, and the remainder in the 
hydrogen electrode. As the set dries out, a larger percentage of the 
electrolyte remains in the nickel electrode. 

The chart on the right describes the percentage of the pore structure that 
is filled with electrolyte. When the cell set is completely filled, then all 
three components are of course filled. As the set of components begins 
to dryout, it is seen that the separator material preferentially loses 
electrolyte due to the separator’s high percentage of large pores 
compared to the other two components. The nickel electrode and the 
hydrogen electrode do not stay fully saturated but do share a small 
amount of their liquid content with the separator. 

Later, it will be shown that a cell displays dryout symptoms when the pore 
structure of the separator is only about 30% filled with electrolyte. At this 
point, the pore structure of the nickel electrode is between 90% and 95% 
filled. 
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2 


CHECKING THE MODEL 


• SET MODEL FOR EITHER FULL CHARGE OR FULL 
DISCHARGE 

• IN THE CHARGED MODE, THE PERCENTAGE OF GAMMA 
PHASE IS INSERTED 

• ALL KNOWN FACTORS - DIMENSIONS, THICKNESSES, ETC. 
ARE INSERTED 

• ASSUMPTIONS RELATED TO COMPRESSION AND 
EXPANSION ARE INSERTED 

• RESULTS ARE NOTED 

- Partitioning of electrolyte 

• RESULTS ARE COMPARED WITH RESULTS FROM DPA 
STUDIES 


Electronics Technology Center tea THE AEROSPACE 
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The model is best set up for the cell to be in the fully charged state or the 
fully discharged state. This permits the proper value for the density to be 
used for the active material. If the cell is set for the fully charged 
condition, one must insert a value for the percentage that is in the gamma 
phase. As the cell approaches a dryout condition, selecting the proper 
density for the active material does make a difference in the amount of 
electrolyte that is calculated to be in the separator. 

Other factors that are known with a higher degree of certainty are then 
added - dimensions, thicknesses, etc. 

Assumptions relating to how much the electrodes are expected to expand 
over the life of the test and how much each of the components are 
compressed due to being compressed between the end plates by the 
Belleville washer are next made. 

The results are noted in terms of partitioning and percentage of the 
separator pores that are filled with electrolyte. 

During the validation process, the results from the model were compared 
with the results from the DPA studies. 
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EARLY ATTEMPTS TO VALIDATE MODEL 


• CELL MODIFIED AND WATER INCREMENTALLY REMOVED* 

• PERFORMANCE MONITORED DURING DISCHARGE AT 
DIFFERENT RATES 

• RISE IN INTERNAL IMPEDANCE WITH REDUCTION IN 
ELECTROLYTE VOLUME WAS NOTED 

• AT THE POINT WHERE THE IMPEDANCE DOUBLED, THE 
CELL WAS DISASSEMBLED AND THE PARTITIONING WAS 
DETERMINED VIA TITRATION OF COMPONENTS 

- THE MODEL PREDICTED A DRYER SEPARATOR THAN THE 
RESULTS OF THE TITRATION ANALYSIS 

‘DETAILS IN 1999 ANNUAL BATTERY CONFERENCE PAPER 
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Early attempts to validate the model consisted of developing a dryout test 
in which incremental amounts of water were evaporated from a known 
new cell (no corrosion or expansion) that had been filled with a known 
amount of electrolyte during its activation. As electrolyte volume was 
removed from the cell, the performance was monitored at several 
discharge rates over the course of the full discharge. The rise in internal 
impedance as the electrolyte volume was reduced was noted. Evidence 
for diffusion limiting currents were watched for, and was occasionally 
seen. When the internal impedance had doubled, the cell was judged to 
have dried out. The cell was disassembled and selected components 
were analysed for their electrolyte content. The results suggested that the 
partitioning of the electrolyte was uniform in the six locations within the 
cell that were examined. The initial versiop of the model predicted a 
dryer separator than was suggested by the chemical analyses. 

Since I was dissatisfied with the agreement between the two sets of 
answers, attempts were made to bring the results and the prediction 
closer together. When this process was completed, the model very 
closely matched the post test analyses. More realistic numbers for the 
compression of cell components had to be determined. Out-of-cell 
thicknesses of expanded electrode material had to be converted to in- 
stack thicknesses. The assumption that the nickel electrodes were 100% 
full gave way to the more valid value of between 90 and 95%. Since the 
cell was taken apart in the discharged state, the density of the beta nickel 
hydroxide had to be used. The remainder of the paper will outline the 
steps in this process. 
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STEPS IN UPGRADING MODEL 


• VARIABLES SET FOR THE DISCHARGED CONDITION 

- Condition when the cells were opened 

• PERCENT EXPANSION IS ADJUSTED TO BEST FIT OF 
ELECTROLYTE PARTITIONING BASED ON DPA 

• ASSUMPTION OF 100% FILL OF POSITIVE PLATES IS 
ADJUSTED TO ACTUAL PLATE EXPANSION 

• PERCENT FILL OF POSITIVE PLATE IS COMPARED TO 
PREVIOUSLY GENERATED CAPILLARY PRESSURE STUDIES 

• THE MODEL RESULTS ARE NOW IN VERY CLOSE 
AGREEMENT WITH DPA RESULTS 
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The chart on the next page will be used to illustrate how the available 
data related to cell dryout was used to upgrade and validate our original 
version of the static model for nickel-hydrogen cells. 
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SOURCES OF DATA FOR THE STUDY 


• RESULTS OF SPECIALLY CONDUCTED DRYOUT TESTS 
WHERE THE ELECTROLYTE VOLUME WAS SEQUENTIALLY 
REDUCED AS THE PERFORMANCE WAS MONITORED 

• POST TEST ANALYSES WHERE THE PARTITIONING OF 
ELECTROLYTE BETWEEN THE SEPARATOR AND THE 
NICKEL ELECTRODE WAS DETERMINED 

• POST TEST ANALYSES WHERE THE AMOUNTS OF 
EXPANSION OF THE NICKEL ELECTRODES WERE 
MEASURED WITHOUT ANY COMPRESSIVE FORCES 

• RESULTS OF COMPRESSING STACKS OF TYPICAL CELL 
COMPONENTS 


Electronics Technology Center laga THE AEROSPACE 

Energy Technology Department IKBlCORPORATION 


We had several sources of information that we were able to use in our 
efforts to upgrade and validate our static model for nickel-hydrogen cells. 
Different pieces of the inputs that were available to us'were complete in 
differing degrees, but they were all helpful. Information on the dryout or 
partitioning aspects of six different cells were used in this study. 

Common to all cases was an evaluation of the cell performance at the 
time the cell was disassembled or the dryout test was completed. The 
three dryout tests were ended at the point where there was a significant 
increase in internal impedance. Full cell capacity was always available, 
but at a reduced terminal voltage. In all cases, information was available 
on the initial amount of electrolyte that was added to the cell. Three tests 
were carried out where the water portion of the electrolyte was 
incrementally removed using our dryout technique. Three data sets were 
available from cycling tests that were stopped and the partitioning of the 
electrolyte between components was determined or estimated. The cells 
covered several different manufacturers, cycling conditions, electrolyte 
concentrations, and nickel electrode types. The compressive 
characteristics of representative cell components from disassembled cells 
were also incorporated into this study. In all cases, it was assumed that 
the results of post test analysis was correct and the results of the static 
model was required to give the same results. All six cells were treated 
equally and the model was expected to be equally accurate when applied 
to each example. 
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COMPONENTS STUDIED 


• GAS SCREENS 

• HYDROGEN ELECTRODES 

• ZIRCAR SEPARATORS 

• NYLON SEPARATORS* 

• UNIMPREGNATED PLAQUE 

• NICKEL ELECTRODES WITH DIFFERENT AMOUNT OF 
EXPANSION - 10%, 20%, AND 30% 

• BELLEVILLE WASHERS* 

• POLYSULFONE CORE PIECE* 


*BASED ON AVAILABLE LITERATURE INFORMATION 


Electronics Technology Center 
Energy Technology Department 


The above listed items are typical of nickel-hydrogen cells. Over the 
years, our laboratory has disassembled many different kinds of IPV 
nickel-hydrogen cells. We were able to select representative 
components that were available to us from these DPA studies. 
Information on different Belleville washers and the tensile characteristics 
of the material used for the core piece were obtained from manufacturer 
data sheets. Earlier SAFT cell designs used nylon-based separator felts 
that had been reported upon at previous NASA Battery Workshop 
meetings.* 


Denning, C. and Jamin, T. “Polyamide Separator Behavior in NiH 2 Cells”, 
Proceedings of the 1994 NASA Aerospace Battery Workshop 
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RESULTS OF COMPRESSION STUDIES 



0 20 40 60 80 100 120 140 160 

Stress (PSI) 


■ Zircar Type 1 screen • Type 2 screen 

▲ Old type plate ▼ New type plate X Expanded 
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The compression studies consisted of stacking eight elements of the item 
under investigation and applying pressures using an Instron machine. 

The amount of displacement was monitored along with the compressive 
forces. From the plot of displacement on a per item basis, it is seen that 
all of the components had the same general shape as they were 
compressed. Some components were more resistant to the compressive 
forces, but they all displayed a tendency to be squeezed to a thinner 
dimension as pressure was applied to them. On this figure, the 
compressibility of three different nickel electrodes are plotted - one taken 
from a cell that was almost new (about 10% expansion relative to the 
original plaque thickness due to the activation step), one that had about a 
20% increase in its electrode thickness, and another had a 30% increase 
in its out-of-cell thickness. These thickness variations resulted from 
different amounts of cycling stress. 

The information in this figure is needed to help estimate the thickness of 
these components while they are still inside the cell under the 
compressive forces of the core piece and the end plates. 

To be noted in the figure is the fact that when the cell is early in life where 
not too much expansion has taken place, the zircar separator is less 
resistive to compression than the nickel electrode. However, after 
extended cycling,. the zircar separator is less compressible than the 
electrode that had expanded about 30%. It is not entirely trivial to 
estimate the thicknesses of all the components while they are still in the 
cell following a life cycle test. 
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ESTIMATING IN-STACK THICKNESSES 


• DETERMINE COMPRESSION CHARACTERISTICS OF ALL 
COMPONENTS INCLUDING THE BELLEVILLE WASHER 

• NICKEL ELECTRODES WITH DIFFERENT AMOUNTS OF 
EXPANSION MUST BE TESTED 

• DETERMINE THE ELONGATION OF THE CORE PIECE 

• BY TRIAL AND ERROR, DETERMINE WHERE THE SUM OF 
INCREASES EQUALS THE SUM OF THE DECREASES 

• THESE DIMENSIONS ARE THEN TESTED IN THE MODEL 


Electronics Technology Center Bfca THE AEROSPACE 

Energy Technology Department Is9 C O R P O R AT I O N 


When cells are built, the plate pack of electrodes, separators, and gas 
screens are compressed between endplates and one or more Belleville 
washers. After the proper degree of compressive forces have been 
applied to the plate pack, the end nut is attached to the central core piece 
of the cell assembly. At this point, all of the components are squeezed 
together. Depending on the compressive characteristics of each 
component, it will be somewhat thinner than it was in the uncompressed 
state. The Belleville washer will be compressed and the polysulfone core 
piece will be stretched. 

As cycling proceeds, the nickel electrodes will expand depending on the 
amount of stress placed on them by the cycling conditions. This will 
result in more compressive forces on the remainder of the components. 
At the end of the test, the components are set free of the compressive 
forces and then measured for the amount of expansion they have 
experienced. The real question is what was their thickness when they 
were still in the cell. It is under the compressed condition in the plate 
pack that the electrolyte is distributing itself between the different 
components. 
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SUMMARY OF RESULTS 


Cell Identification 

“Cell A - 

“CeUB" 

Cell'C" 

"CeiTD - 

Cell E 

CellT - 

Cell Dryout Test 

Yes 

Yes 

No 

No 

Yes 

Yes 

Fill Concentration -% 

38 

31 

31 

31 

31 

26 

Partitioning - by analysis 







In Separator % 

20 

20 

37 

42 

N/A 

N/A 

In Nickel Elect % 

77 

75 

60 

54 

N/A 

N7A 

Partitioning - from model 







h Separator % 

20 

15 

37 

42 

22 

6 

In Nickel Elect. % 

75 

73 

59 

54 

74 

”90 1 

In Hydrogen Elect % 

5 

6 

4 

4 

5 

4 

Fill % - Based on Model 







Overall Cell % 

80 

74 

77 

82 

63 

71 

Separator % 

45 

30 

56 

65 

27 

13 

Plate Expansion 







From DPA-% 

12 

6 

20 

12 

12 

N/A 

Used in Model-% 

6 

0 

8 

5.5 

6 

0 

Corrosion of Ni -% 

0 

0 

6.3 

5.5 

0 

0 








% Vol. Loss at Dryout 

N/A 

25 

N/A 

N/A 

20/32 

28 
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The information that was available to us was placed in the table above. 
The initial conditions and fill amounts were all known and the out-of-cell 
plate thicknesses were measured at the end of the test: Cells A through 
D had measurements of the partitioning between components. At this 
point, the static model for each of these cells was set for full discharge, 
the separator was assumed to be 10% compressed, the amount of 
corrosion was inserted, the amount of plate expansion as measured out- 
of-cell was inserted, and the model was solved. 

If the results did not match the results of the partitioning that were 
obtained by chemical analysis, the amount of expansion of the nickel 
electrodes was adjusted so that a match occurred. This process worked 
for three of the four cells. Cell B would have required a net compression 
to bring the model results in line with the analysis data. With the 
exception of cell C, an adjustment of about 6% in the as-measured 
expansion was necessary to bring the results of the static model in line 
with the results from the chemical analysis. 

Following this, the assumption used in the model which stated that the 
nickel electrode was always 100% filled was changed to be in keeping 
with the more accurate results of the earlier NASA pore size studies. By 
assuming that only 95% of the nickel electrode pores were filled with 
electrolyte rather than 100%, it was no longer necessary to apply any 
adjustments to the electrode thickness. Irrespective of what the model 
predicts, the evidence from the different dryout and partitioning studies 
demonstrate that when the separator is only about 30% filled with 
electrolyte, significant performance degradation results due to dryout. 
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SUMMARY 


• THE LITERATURE HAS DOCUMENTED A FEW CASES OF 
CELL DRYOUT 

• SPECIAL DRYOUT TESTS AND THE AVAILABLE 
INFORMATION RELATED TO PLATE EXPANSION AND 
SUBSEQUENT PARTITIONING OF ELECTROLYTE SUGGEST: 

- Dryout results in significant increases in resistance before it 
results in diffusion-limited currents 

- Double layer Zircar separators are almost impossible to dry out 

- Single layer Zircar designs should not be rejected out of hand 

- The original static model predicted a dryer separator than that 
found in DPA studies 

- Upgraded model useful for evaluating cell designs and projected 
cycle life based on known degradation modes 


Electronics Technology Center flirifl THE AEROSPACE 

Energy Technology Department ifislCORPORATION 


There have been documented and undocumented cases of cell dryout 
during life cycle testing. There have also been documented and 
undocumented cases of excessive amounts of expansion of nickel 
electrodes leading to breakage of the core piece and/or failure of the 
threads between the core piece and the core nut. Knowledge of the 
impact of these phenomenon is helpful in finding ways to avoid these 
situations. The static cell model which has now been updated and 
validatied coupled with estimates for corrosion and expansion rates as 
affected by cycling conditions, can be used to assess either a suggested 
cell design or a suggested set of cycling conditions and recharge 
protocols. 

Some cell designs use two layers of separator material, whereas others 
use only one. The weight savings resulting from using only one layer are 
significant. The model coupled with the understanding of the changes 
taking place during cycling can be used to evaluate the suitability of using 
the lighter-weight single layer cell design. 

As dryout is approached, a small error in the accuracy of the model 
resulted in significant changes in the suggested percentage of the pores 
of the separator that were filled with electrolyte. For this reason, it was 
necessary to be as accurate as possible in terms of the assumptions that 
were used in the model. Further, these assumptions had to be checked 
against the actual data that were available following different kinds of 
tests that addressed the issues of cell performance and the partitioning of 
electrolyte between components. We now feel very comfortable with the 
accuracy and utility of our nickel-hydrogen model. 
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CRANE Cycles Ongoing 
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Depth of Discharge (%) 





(%) Ap«dB3 HV aiqBSfi 
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Cycle Life (%) 




1999 NASA Aerospace Battery Workshop 


- 560 - 


Nickel-Hydrogen Session II 




( J HAV> 

SuiuiBiusy 

3 

Q 

30 

d 

B 

K 

S: 

r 

m 

B 

d 

4\ 

0S 

(N 

30 

(N 

(HV) 
Xjidb(1b3 |bjox 

N 

3 

1 

r 

o 

r 

d 

2 i 

2 

r- 

n 

n 

V" 

(HV) 

aSesfi ejpQ 

30 

d 

rt 

t" 

00 

ft 

ro 

d 

3 ( 

d ■ 

rr ' 

00 

d 

d 

(HV) XjpBdBO loa 
* aoa 3a v 

■o 

O 

d ! 
d 1 

»r, 

© 1 
d i 
*n i 

r- i 

00 - 
in • 

rn 

fn 

0 

(%) XjpBdBO 
aiqespi SmmBUJOfl 

30 

3S ! 
OS i 

Os 

- 

rd i 

d i 

O i 

00 1 
00 < 
OS i 

ro 

d 

■0 

0 

d 

■0 

(%) P»pu»dx3 
3jn wro 

o 

d 

2 

c* 

3 

5i 

so < 

m 

d 

OS 

1/3 

d 

cs 

JO Pd 

s 

d 

<N 

d 

<N 

«r, 

<N - 

0S 

d 

(N 

rM 

d 

(mg |enpiAipui) 
^Hiqede;) 3jn 

43,188 

48,948 

r<) « 

-r 

cs 

d 

-r 

41,573 

r- 

os 

0 

00 

m 

V) 

» n 
00 

d 

T 

(%> 
aoa 8 av 

o 

d 

so 

<N 

tTi 

1 

< 

d < 

i n ' 

3 

d 

0 

0 

d 

so 

(%) 

aoa “la 

o 

d 

so 

rd 

> n 

i n ' 

O < 
ri \ 
O ' 

30 

d 

O 

00 

d 

*n 

aSjBipsfQ 

o - 
— 

rk 

o - 
o < 
© - 
<N t 

3 < 
3 < 
rs 

<N - 

3 

3 

-r 

M 

3 

3 

O 

N 

3 

0 

<N 

uoissiiv 

dSjBqosjQ 


o < 
o - 

d « 

3 < 
3 < 

-r 
3 < 

3 < 
3 < 
TT ■ 
3" < 

3 

3 

-r 

3 

§ 


- 561 - 


Nickel-Hydrogen Session II 


1999 NASA Aerospace Battery Workshop 


Using the cycle life vs. DOD and capacity fade characteristics, the remaining energy is 
calculated for each energy bin. This process takes into account the capacity 
degradation caused by each energy bin. The requirement is that the remaining energy 
must be sufficient to support the energy bin at an acceptable DOD. 
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In this case the state of charge will rundown during the orbital day. 
This introduces the additional requirement that the degraded capacity 
must support each energy bin starting at less than full charge. 
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The modeling table has been updated to incorporate an “adjusted” energy 
requirement to accommodate the rundown condition. 
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Positive Precharge Characterization 
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Experience base 
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precharge configuration, provide certain advantages and are 
being baselined for many missions. 
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NiH 2 Cell Chemistry 
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Verification of Nickel Precharge 
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Precharge Verification 

Nickel Precharge - Eagle Picher RNH160-3 Cells 
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Nickel Precharge -- Eagle Picher RNH98-1 Cells 
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Strain Gage Effects and the Van der Waals’ Equation 


By 


Dean W. Maurer 


Introduction 


One of the many advantages of the nickel hydrogen battery is that it provides a direct and 
continuous measure of its available capacity. Since hydrogen gas is the active material of the 
negative electrode and its pressure is a direct measure of its quantity (with certain restrictions to be 
discussed later), a simple measurement of the pressure and temperature can be converted to 
capacity available (assuming an offset for the pressure remaining at cut-off). The pressure is 
measured with a strain gage which is bonded to the dome of the pressure vessel and responds to 
the microscopic reversible bulging of the dome with pressure (the cell wall is only 0.024” thick). 

In addition since hydrogen gas obeys Faraday’s law and well known equations of state (gas laws), 
measures of the rate of change of pressure can be used to derive other parameters such as power 
levels or currents during charge and discharge. In attempting to do this for Telstar batteries, it 
became apparent that the strain gage output itself varied with temperature. This paper describes 
methods used to determine the cell free volume and to calibrate the strain gages for temperature 1 
effects. 


The Strain Gage 


The strain gage consists of two stress sensitive resistors wired into the arms of a Wheatstone 
bridge. One of them is rigidly bonded onto the dome (active element) and the other (the dummy) 
is bonded along with the rest of the circuit to a piece of metal from the dome of another cell can 
which in turn is bonded to the dome. The other two resistors of the bridge are adjusted to create an 
output signal in the appropriate range (~10 mv max) when the bridge is activated with the 
specified excitation voltage (typically 5 to 10 volts). This use of the two stress sensitive resistors is 
intended to compensate for temperature effects such as metal expansion. After the cell is built, but 
before electrolyte addition, each strain gage is calibrated at room temperature at 4 or 5 different 
pressures against an external calibrated dial gage. Typically these data fall on a straight line with 
approximately a 0.01 mv/ psi slope and an intercept close to zero and a coefficient of correlation 
greater than 0.9999. It is important to keep in mind that the strain gage measures the gage 
pressure not the absolute pressure. In practice when the cell pressure falls below atmospheric 
pressure the strain gage output remains constant or nearly so. For calculations described here 14.7 
psi must be added to the value provided by the strain gage. Also note that when the cell is in 
space the gage and absolute pressures are the same. This needs to be considered in those 
circumstances. 
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The Equations of State 

The simplest equation of this type is the Ideal Gas Law: 

PV - nRT ( 1 ) 

Where P, V and T are the absolute pressure, volume and absolute temperature. ”n” is the number 
of moles 1 of the gas and R is the gas constant equal to 0.082059 liter-atmospheres/ mole-°K or 
1206.238 cc-psia/ mole-°K. In this form the molecules are considered to be point masses with no 
interaction other than kinetic collisions. It holds very well at low pressures and for simple 
molecules; for hydrogen it is good up to about 200 psia. For higher pressures the more elaborate 
Van der Waals’ equation is required: 


P = 


nRT 

(V -nb) 



VV 


y 

) 


( 2 ) 


Where the variables are the same as before and the constants a and b are the Van der Waals’ 
constants which depend on the gas. For hydrogen, a = 0.2444 liter 2 -atmosphere/ mole 2 or 3592680 
cc 2 -psi/ mole and b = 0.02661 liters/ mole or 26.61 cc/mole. In this equation, the term with “a” is 
an adjustment to the pressure to account for the attraction between molecules due to their close 
proximity (Van der Waals’ forces). This attraction is the result of electric dipoles in the molecules 
caused by the molecular distribution of electrons. For hydrogen this is a small effect; it would be 
larger with a polar molecule such as water or carbon dioxide. The term with “b” is a correction to 
the effective volume to account for the finite volume of the molecules themselves. The two 
equations give similar results below about 200 psia but gradually diverge above that pressure and 
are about 5% different at 800psia. This difference will be the topic of the next section. Equation 
(2) can be rearranged as: 

PV* - n(bP + RT)V 2 + an 2 V-abn* =0 (3) 

Note equation (3) is cubic in both V and n. Fortunately the numbers are such that this equation can 
be solved with conventional techniques and has a real root. In Excel spreadsheets the solution can 
easily be written as a “function macro” in Visual Basic. 

The Van der Waals/ Ideal Gas Ratio or (Pi/P v ) 

Using equations (1) and (2) above the Pj/P v ratio can be calculated for specific values of V and T 
and a series of n’s. A plot of the ratio vs. P v results in a quadratic form. The curve fit values of sets 
of A’s and B’s for each set of particular V’s and T’s (the C value of the quadratic expression is 1 
because at low pressure the two laws converge). By repeating this operation for the useful range of 
values of V and T, a series of A’s and B’s are obtained which can in turn be curve fit by 
multivariable least squares. The net result for V in the range of 900 to 1200 cc’s and T in the range 
of 250 to 3 10° K is: 

Pi /> = 1 - (7.48 1 8 * 10 -9 + 8.7204 *10" 13 *K-2.7031*10 _ " *t)*Pv 2 
/Pv ' 1 (4) 

- (6.4571 * 1 0 -5 - 8.393 * 10~ 10 * K -5.707 *10' 8 *f)*Pv 


1 A mole of a substance is a specific amount which has a weight equal to its molecular weight and is 

6.023x1 0 2 , molecules or Avogadro’s number of molecules. 
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In this equation the volume is in cc’s and the temperature is in degrees Kelvin. This equation can 
also be written as function macro. 

Another technique is to rearrange equation (2) above and to divide through by V, which gives: 

(F - nb^Pv + a(n/V) 2 ) nRT 

V v K 

The term to the right of the first equal sign is seen to be the ideal gas expression and equals the 
ideal gas pressure Pi. Equation (5) is more generally applicable than equation (4), however, in the 
applicable range of (4), it is more useful if it is inconvenient to use “n” whereas (5) is more useful 
(simpler) when “n” is available. 

Faraday’s Law & Graphical Techniques 

In simple terms Faraday’s law states that when 96,496 Coulombs or amp-seconds of charge are 
passed through an electrolysis cell, one equivalent of matter (Avogadro’s number of electrons or 
charges) is reacted at each electrode. The 96,496 Coulombs is called “The Faraday” (F) and is 
usually rounded off to 96,500 Coulombs. This concept, using electroplating of silver under 
specific conditions, defines the “ampere”. For the hydrogen electrode as we use it, there are 
essentially no parasitic reactions on discharge 2 so the relationship between charge passed and the 
quantity of hydrogen is considered to be exact. In units of amp-hours (AH) the Faraday is 26.3 
AH/ equivalent 3 . Since for the hydrogen reaction there are two equivalents per molecule, this 

results in 53.6 AH/ mole of hydrogen. Thus, 

n = At/53.6 (6) 

where t is time in hours and A is the current in amperes. Combining equation (1) with (6), gives: 


Pi = 


AtRT 

53.6V 


Thus if Pi is plotted vs. t and the slope determined then, 


(7) 


dPi 


= slope = 


ART 


( 8 ) 


dt 53.6V 

from which one of the constant terms can be determined if the others are known. In addition, if T 
is not constant Pi/T can be plotted and A or V determined. If the discharge is also at constant 
power, W/E= A (where E is the battery voltage) can be substituted into the equation so that from a 
plot of EPi/T vs. t, W can be determined from the slope. 


The Cell Volume 


The free volume of the cell can be determined by a variety of techniques but the selection is 
limited after the cell is sealed. One method is to use the well-controlled, constant temperature 
capacity tests, available in Excel, made at the supplier as part of the cell lot acceptance tests 
(ATP). Starting with the 20°C data since the strain gages (SG’s) were calibrated at room 
temperature, the pressure (after converting it to absolute pressure) is plotted vs. time starting from 
about 200 psi down to the 1 .0 volt cut-off. From the slope, the volume can be estimated (equation 
(8)). This value has some errors associated with it which will be discussed bellow. 


2 Actually there is some self-discharge but this is negligible at the normal discharge rates and temperatures 
used. The charge reaction does have parasitic elements that are significant, namely oxygen evolution, self- 
discharge and the direct reaction of hydrogen with the charged nickel species on the nickel electrode. 

J That is per electron transfer 
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Next from the current and discharge time a data column is generated (Table 1) providing the 
capacity removed at each data point. Using the volume estimate and the ideal gas law the 
44 hydrogen equivalent capacity” (HEC) remaining at the cut-off (or any other point) can be 
estimated from the absolute pressure reading at that point. From this value a column of U HEC 
remaining” at each point in the discharge can be generated. Using a macro generated to determine 

Table 1 A Sample Section of Spreadsheet Data 


Discharge 

Time 

Capacity 

Out 

Average 

Temp 

HEC 

Remaining 

Supplier 
SG Pres. 

Apparent 

Volume 

Back Calc. 
SG mv 

VdW 

Pressure 

0 

0 

19.3 

121.8 

816.35 

1003.76 

8.53 

817.07 

0.167 

12.52094 

19.9 

109.2 

732.01 

1000.23 

7.63 

729.92 

1.167 

87.52499 

18.5 

34.2 

213.43 

996.54 

2.13 

211.96 

1.333 

99.97621 

19.4 

21.8 

126.08 

1027.27 

1.20 

129.53 

1.417 

106.2756 

18.5 

15.5 

87.36 

1002.94 

0.79 

87.36 


AVE T = 

19.1 

15.5 

= CUT OFF 
HYDROGEN CAP. 




the volume from cubic solution of equation (3), the calculated apparent volume at each step of the 
discharge is determined. An example of a plot of the data for a cell is shown in Fig 1. A trendline 


Fig 1, Apparent Volume of S/N1295 @ 20 degC 

y = -4.75x + 1002.95 


LU 

a: 
< 
CL 
Q. 





Discharge Time (hrs) 


is generated and the intercept is used as the best representation of the free volume of the cell. This 
volume will be different from the one used earlier to determine the HEC at cut-off. These steps are 
repeated using the new apparent volume estimate to recalculate HEC at cut-off and recalculate the 
column of “HEC remaining” at each step. The apparent volume and Figure 1 is redone and a new 
volume determined. These steps are repeated as many times as necessary for the values to 
converge; usually about five times. This effectively averages a number of data points but note that 
the trendline has a slope. The possible causes of this slope will be discussed in the next section. 

This process was repeated for the other ATP temperatures with a surprising result, namely a 
systematic variation of the apparent volume with temperature as shown in Fig. 2. Since changes in 
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the cell free volume could be caused by volume expansion of the cell components, differences in 
density of the active material, etc., changes of this magnitude are physically unreasonable, other 
possibilities were explored, namely the strain gage calibration. 



Strain Gage Temperature Calibration 

Returning to the ATP data in Table 1, the column generated previously for the ”HEC remaining” 
is used to calculate the absolute pressure with a macro based on the Van der Waals equation (2). 
The volume previously determined in the 20°C test, is used for this calculation since the strain 
gages were originally calibrated at that temperature. This absolute pressure is converted to gage 
pressure and plotted vs. the strain gage voltages derived from the gage pressures using the original 
supplier calibrations. A trendline then gives the new calibration values. Figures 3 & 4 show these 


Fig. 3 SG Temperature Effect 



T em peratu re 


y = -0.1765x ♦ 145.7 
R 2 = 0 9994 


y = -0 1 528x + 142 9 
R 2 = 0 9993 
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effects on both the slope and the intercept. Note that the slope changes by 8% in 40 degrees or 0.2% per 
degree. The intercept is non linear and varies by 15 psi in 40 degrees or the equivalent of about 2.4 AH for 
cells with lOOOcc of volume. Using the two equations from the trendline fits a revised strain gage 
calibration equation can be generated in which the slope is itself a linear equation in T and the intercept is a 
quadratic equation in T. In cases where it is inconvenient to use this approach many times the apparent 
volume equation from Figure 2 can be used instead. It actually embodies the same effect even though it is 
physically incorrect. The cause of the slope in the apparent volume data in Figure 1 can now be explained 
with the additional fact that thermocouple measurements in the strain gage region of the cell during high 
rate discharge show a significant temperature rise. The higher gage temperature is reflected in a lower 
apparent pressure reading, which causes the slope of the P vs. t plot to be greater. From equation (8) it can 
be seen that a higher slope will give a smaller volume. 


An example of the use of these data to predict the discharge current is shown in Figure 5 which is a plot of 
the derived ideal pressure, equation (4), vs time for the 10°C test. From the slope, the current is calculated 
using equation (8) which gives a value of 74.93 A or within 0.09% of the actual measured value. 


Fig 5 Ideal Slope @ 10 deg C 
S/N 1295 



Slope Current = 74.93 
Measured Current = 75.0 


I 


Discussion 


Manipulating strain gage data is still a work in progress. Several comments can be made: 

1 ) If no allowance is made for the effects described above (supplier calibrations are used as is) 
an error of 5% or more may exist especially at temperature extremes. The capacity estimate at 
the low temperature at the end of charge will be understated (conservative) and the capacity 
remaining near the end of discharge (higher temperature) will be overstated (optimistic). 

2) Cursory review of application of these equations to situations involving more rapid 
temperature changes or to capacity retention tests has revealed large errors; for example in a 
ramped discharge temperature case. 

3) It appears that there may be problems with local internal temperature gradients. That is during 
high rate discharge the cell wall temperature may be quite different from the temperature of 
the bulk of the hydrogen. 

4) The charging process needs to be addressed. 

Conclusion 

The strain gage output is seen in this paper to have many peculiarities and its use to estimate capacity 
available in a battery needs careful attention if accuracy greater than 5% or so is required. 
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